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ABSTRACT 

The origin of the observed time lags, in nearby active galactic nuclei (AGN), between hard 
and soft X-ray photons is investigated using new XMM-Newton data for the narrow-line Seyfcrt I 
galaxy Ark 564 and existing data for 1H 0707-495 and NGC 4051. These AGN have highly 
variable X-ray light curves that contain frequent, high peaks of emission. The averaged light 
curve of the peaks is directly measured from the time series, and it is shown that (i) peaks occur 
at the same time, within the measurement uncertainties, at all X-ray energies, and (ii) there 
exists a substantial tail of excess emission at hard X-ray energies, which is delayed with respect 
to the time of the main peak, and is particularly prominent in Ark 564. Observation (i) rules out 
that the observed lags are caused by Comptonization time delays and disfavors a simple model 
of propagating fluctuations on the accretion disk. Observation (ii) is consistent with time lags 
caused by Compton-scattering reverberation from material a few thousand light-seconds from 
the primary X-ray source. The power spectral density and the frequency-dependent phase lags 
of the peak light curves are consistent with those of the full time series. There is evidence for 
non-stationarity in the Ark 564 time series in both the Fourier and peaks analyses. A sharp 
'negative' lag (variations at hard photon energies lead soft photon energies) observed in Ark 564 
appears to be generated by the shape of the hard-band transfer function and does not arise 
from soft-band reflection of X-rays. These results reinforce the evidence for the existence of 
X-ray reverberation in type I AGN, which requires that these AGN are significantly affected by 
scattering from circumnuclear material a few tens or hundreds of gravitational radii in extent. 

Subject headings: galaxies: active - galaxies: Seyfert - Xrays: galaxies - Accretion, accretion disks 
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1. Introduction 

Active Galactic Nuclei (AGN) in the X-ray 
regime can be highly variable over a wide range 
of timescales. Two common statistical tools for 
analyzing this variability are the power spectral 
density (PSD) of the time series, and the lag spec- 
trum, defined as the time lags between light curves 
obtained from two bands of photon energy and an- 
alyzed as a function of the frequenc y of Fourier 
modes (e.g. iNowak fc Vaughanlll996l ). A positive 
lag is defined as the hard band lagging the soft 
in a particular frequency band. The lag spec- 
trum provides an important constraint on models 
that seek to explain the observed variability (e.g. 



Nowak fc Vaughan 1996; Miller et al. 2010 



lty 
Ja) 



The first clear detection of both positive and 
negative lags in an AGN was in 1H 0707-495 
( Fabian et al.l 120091 ). At low frequencies (v < 
5 x 10~ 4 Hz) the AGN exhibited a positive lag, 
but at higher frequencies the lag was persis- 
tently negative. Similar results have been ob- 
tained in other AGN (e.g. iMiller fc Turnerl 12011 



Emmanoulopoulos et al. 201 it De Marco et al 



2012). The same features are observed in our 



analysis, prese nted here, of new obs ervations of 
Ark 564, which iMcHardv et all (|2007h have previ- 
ously suggested displayed both positive and neg- 
ative lags from fits of Lorentzian components to 
the PSD and cross-spectrum. 

A number of models see k to explain t he la; 
spec trum of 1H 7 07-49 5 
(201 



Fabi an et al.l (200 



and Zoghbi et al. 

lArevalo fc Uttlevl (Ex) 



applied the model of 
invoking propagating 
fluctuations in the accretion disk to explain the 
positive lags at low frequencies, together with 
reverberation of X-rays reflecting from within a 
gravitational radiu^ of a rapidly rotating black 
hole to expl a in the high-frequency negative lags. 



Miller et al.1 (|2010al ) suggested reverberation due 



to scattering from gas 10s to 100s of gravitational 
radii from the central black hole, showing how 
both positive and negati ve lags may arise from 
the same process (see also Miller fc Turnerll201lh . 
Older work on the positive lags seen in galactic 
binary systems and AGN also investigated the role 
of Comptonization up-scattering time delay s (e.g. 
Nowak fc Vaughan|[l99l INowak et al.lll999t) . 



1 Light travels a distance of one gravitational radius (r g = 
GM BH /c 2 ) in 10 s for black hole mass M BH = 2 X 10 6 M Q . 



Analyses of X-ray AGN time series have usu- 
ally been restricted to Fourier space, implicitly 
assuming that AGN time series are statistically 
weakly stationary an d often assuming t hat they 
are Gaussian (but see Uttlev et al. 2005). All the 
above proposed explanations may be considered 
as examples of a 'moving average' time series (e.g. 
Brockwell fc David 120021 ) . in which an underly- 



ing white noise process has been convolved with 
a transfer function (considering only the intrin- 
sic source variations free of measurement noise). 
But the Fourier methods used to date only extract 
limited information on that transfer function. It is 
possible to obtain a good estimate of the PSD, and 
of the lag spectrum between two energy bands, but 
it is not possible to extract the phase information 
of the transfer function in the time series for any 
individual energy band. This information is vital 
for reconstructing the transfer functions, and for 
testing the various models. The aim of this paper 
is to approach the problem in the time domain, 
and directly extract the transfer functions from 
the individual time series by examining the most 
prominent flares in the light curves. 

In this study, we examine three highly X-ray 
variable AGN: the nearby narrow-line Seyfert I 
galaxies Ark 564 (redshift z = 0.0247), 1H 0707- 
495 {z = 0.0411) and NGC 4051 (z = 0.0023). 
NGC 4051 has a black hole mass A/ B h = 1-73! 



.55, 
.52' 



10 M de termined by optical reverberation mea- 
surement ( Dennev et al.l 2009). The other two 
AGN do not have optical reverberation mea- 
surements published at the time of writing. 
Botte et al.l (|2004|) have estimated the mass of 
Ark 564 as M BH = 2.61 ± 0.26 x 10 6 M o using 
the estimated relation between the radius of the 



broad - li ne region and l umino sity of iKaspi et al 



(|2000|). IZhou fc Wand (|2005f) have also used the 



Kaspi et al.1 relation to estimate the mass of 
1H 707-495 as M R H ^ 2.3 x 10 6 M , but we note 
that iLeighlvl ((200J) infers that the mass would be 



M B h ~ 1O 7 M if 1H 0707-495 were radiating at 
its Eddington luminosity. The ultraviolet spectra 
of 1H 0707-495 have bee n ext ensively mode led 
bv lLeighlv fc Moord (|2004l ) and |LeighlyJ (|2004l) as 
being dominated by an outflowing, optically-thick 
wind from the accretion disk. The data used in 
this pap er for 1H 0707-495 is the same time series 
used by IMiller et al.l (|2010a| ). For NGC 4051 w e 
use the Suzaku data used by IMiller et al. ( 2010b ). 
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The Ark 564 time series is new data from XMM- 
Newton. 

2. Observations and Data Reduction 
2.1. Ark 564 

For Ark 564 we obtained a set of eight XMM- 
Newton pointed observations, made between May 
and July 2011. The OBSID, the starting date, 
and the exposure duration of each pointing are re- 
ported in the first three columns of Table [T] Given 
the large energy bandpass and high effective area, 
we focused our analysis on the European Photon 
Imaging Camera (EPIC) pn data. All the EPIC 
pn observations were performed in Small Window 
mode and used the thin optical blocking filter. 
Data were processed using HEAsoft v. 6. 12 and 
SAS v. 11. 0.0. 

The raw pn event files were filtered retaining 
only best-quality (FLAG=0) singles and doubles 
(PATTERN < 4) events. Source events were ex- 
tracted from circular regions with a 36" radius, 
while background events were extracted from two 
boxes of 61" x 100" and 51" x 51" size. Background 
light curves were then extracted for single events 
only, in the 10 — 12 keV band, for each observa- 
tion. After an inspection of these light curves, 
strong flaring background time intervals were re- 
moved using a threshold of 0.1 count s _1 . The 
net exposure times after the flaring background re- 
moval is reported in the fourth column of Table Q] 
These values also account for the Small Window 
mode live time of ~ 71%. 

The SAS task epatplot was then used to esti- 
mate the fraction of pile-up affecting the data in 
the 1 — 9 keV band. We found non-negligible de- 
viations (of the order of 3 — 5%) of the observed 
pattern distribution from the predicted one in half 
of the observations. For the affected observations, 
the cores of the source distribution on the im- 
age were excised, using a radius dependent on the 
amount of pile-up in the data, until we were able to 
recover the predicted pattern distribution within 
a 1 — 2% uncertainty. For the observations where 
we needed to excise the core (i.e. where we used 
an annular source extraction region), the source 
count was then corrected by a factor that repre- 
sented the point-spread-function correction from 
the excised annular region to the count rate one 
would have obtained from a circular source ex- 



traction region with no excision. This correction 
factor ranged from 0.45 — 0.75 and was constant 
within each individual observation (Columns 6-8 
of TableQ] show count rates corrected by this fac- 
tor). Finally, energy dependent exposure maps 
were generated through the SAS task eexpmap 
and the background count rates were corrected for 
the slightly lower effective exposure found relative 
to the source region, and also scaled appropriately 
for each source extraction cell, before the back- 
ground rates were subtracted from the source light 
curves. 
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Fig. 1. — Ark 564 light curves in the soft (0.4-1 
keV, left) and hard (4-7.5 keV, right) bands for the 
eight observations (earliest at top). Light curves 
have been binned to 128s per time bin, corrected 
for excision, and scaled relative to the average 
count rates of 25.7 counts/s and 0.60 counts/s for 
the soft and hard bands respectively. Scaled 68% 
confidence regions are too small to be marked on 
the plots, but have mean values of 0.026 and 0.17 
for the soft and hard bands respectively. 
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Table 1: Summary of 2011 XMM-Newton observations of Ark 564. 



OBSID 


start date 


duration 


exposure 


excision 


0.4-1 keV rate 


4-7.5 keV rate 


0.4-10 keV rate 






/ks 


/ks 


n 


As- 1 






0670130201 


2011 May 24 


59.5 


38.6 


8 


33.51 ±0.04 


0.752 ±0.007 


48.79 ±0.04 


0670130301 


2011 May 30 


55.9 


34.3 





21.67 ±0.03 


0.479 ± 0.004 


31.33 ±0.03 


0670130401 


2011 Jun 06 


63.6 


31.1 





21.56 ±0.03 


0.562 ± 0.004 


31.84 ±0.03 


0670130501 


2011 Jun 11 


67.3 


40.6 


5 


26.84 ±0.03 


0.643 ±0.005 


39.20 ±0.05 


0670130601 


2011 Jun 17 


60.9 


33.5 





22.89 ±0.03 


0.565 ± 0.004 


33.60 ±0.03 


0670130701 


2011 Jun 25 


64.4 


29.1 





14.17 ±0.02 


0.356 ±0.004 


20.69 ±0.03 


0670130801 


2011 Jun 29 


58.2 


40.5 


4 


23.47 ±0.03 


0.570 ± 0.004 


34.03 ± 0.04 


0670130901 


2011 Jul 01 


55.9 


38.7 


8 


30.74 ±0.04 


0.678 ±0.007 


44.36 ± 0.04 



Note. — Columns show observation ID, start date, duration of the observation, net exposure time after live time correction 
and flaring background removal, excision radius used to correct for pile-up and average count rates after correction for excision, 
in the 0.4-1 keV, 4-7.5 keV, and 0.4-10 keV bands. 



Light curves after correction and scaling are 
shown in Fig. [I] The count rates have been cor- 
rected for the excision performed, and the light 
curves in each band have been scaled relative 
to the overall average count rate. The average 
corrected count rate for the 0.4-1 keV band is 
25.7 s -1 , and for the 4-7.5 keV band is 0.60 s -1 , 
with typical 68% confidence regions of 0.026 and 
0.17 respectively. 

2.2. 1H 0707-495 

XMM-Newton observed 1H 0707-495 during 
four epochs, 2000 Oct, 2002 Oct, 2007 Jan-Feb 
and 2008 Jan. In this paper we present analy- 
sis of only the larg est dataset, from 2008, as in 



Zoghbi et al. (201 



irg' 
10) 



In this case, the EPIC observations were made 
using the medium filter and with prime large win- 
dow mode during 2008. These data p roducts were 



based upon the original analysis of I Miller et al 



(|2010aD . processed using SAS v8. 0.0 and HEA- 
SOFT v6.8. The pattern selections used were the 
same as those for Ark 564. The background levels 
were high d uring some of the 1 H 0707-495 obser- 
vations. See lMiller et al.l (|2010al ) for details of the 
filtering of background flares, and the selection of 
source and background extraction cells. 

2.3. NGC 4051 

Suzaku observed NGC 4051 during 2005 Nov 



10-13 and 2008 Nov 6-12 and 23-25. Here we ana- 
lyze only the relatively unobs cured 2008 observa- 
tions fseelMiller et al.l l2010bf). The Suzaku X-rav 



Imaging Spectrometer (XIS. iKovama et al J 12007) 
is an ensemble of CCDs that yield useful data over 
~ 0.6 — lOkeV, with energy resolution ~ 150 eV 
at 6keV. Here we use data from the two front- 
illuminated units, XIS and XIS 3. 

In brief, the data were reduced using v6.4.1 of 
HEASOFT. They were screened to exclude peri- 
ods during and within 500 seconds of the South 
Atlantic Anomaly, to have an Earth elevation an- 
gle less than 10° and to have a cut-off rigidity 
> 6GeV. We selected event grades 0,2,3,4 and 
6 and cleaned out the hot and flickering pixels. 
Source products were extracted from circular re- 
gions of 2.9' radius with background spectra from 
a region of the same si ze, offset from the so urce. 
For further detail see Turner et al.l ([2009?) and 
Miller et al.l (|2010bl ). 



3. PSD and lag spectrum analysis 

In this paper, the time series have been anal- 
ysed in four bands of photon energy: 0.4-1 keV, 
1-2 keV, 2-4 keV, 4-7.5 keV, similar to those used 
in earlier papers. To test models of reverberation, 
it would be good to analyze separate line emission 
and continuum time series, as in optical reverbera- 
tion studies. However, in current X-ray data, pho- 
ton statistics are limited and it is not possible to 
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Fig. 2. — PSDs of Ark 564 (a) (left panel) in 0.4-1 keV (lower circles) and 1-2 keV (upper, squares) bands, 
sampled at frequency resolution A log v = 0.05 (b, right panel) in 0.4-1 keV (lower circles), 2-4 keV (squares) 
and 4-7.5 keV (triangles) bands, sampled at frequency resolution Alogi' = 0.15. Vertical bars indicate 68% 
confidence region uncertainties, horizontal bars indicate the range of frequencies incorporated in each point. 



cleanly separate line emission from continuum: it 
is necessary to work in broad bands of photon en- 
ergy so that variations on short timescales are not 
completely dominated by photon shot noise. 
The PS D of 1H 070 7 -495 ha s previously been 
d2010ah . 



shown in Miller et al 



The PSD of 



NGC 4051 f or this dataset has been shown by 
IMiller et al.l (l2010bl) as w ell as in previous work 
bv iMcHardv et al.l (|2004l). Ark 564 has p revi- 



ously been studied by lArevalo et al.l ((2006) and 
IMcHardv etaD <l2007h . Fig.0 shows the PSD of 
the new observations of Ark 564 in each band, 
ob tained using the ma ximum-likelihood method 
of IMiller et all (l2010al lbh which corrects for the 
gaps in the time series caused by breaks between 
observations and the filtering out of bad data dur- 
ing observations, and also for the presence of shot 
noise. The PSDs are normalized to the mean 
of the full series. The method provides accurate 
estimates of the PSD uncertainties and of the 
covariance between frequency channels, although 
the latter arc difficult to convey in simple dia- 
grams. The two softest bands have high count 
rates and their PSDs arc shown at high frequency 
resolution in Fig.[2][a). At high frequencies the fre- 
quency channels arc largely independent and have 
widths A log 10 v = 0.05. The maximum frequency 
shown is the Nyquist frequency corresponding to 
the chosen time series sampling of 128 s. At low 
frequencies the finite duration of the individual 
observations causes nearby Fourier modes to be- 



come correlated, so the frequency bin widths have 
been increased to preserve approximately the sta- 
tistical independence of the plotted points. The 
PSDs of the two harder bands are shown at lower 
frequency resolution in Fig.[2jb), because of the 
larger statistical uncertainties. We also show the 
0.4 — IkeV band at this resolution, for reference. 

We note in the Ark 564 PSDs the features previ- 
ously reported by lArevalo et al.l : IMcHardv et all 
namely the very steep fall-off to high frequency 
and the complex shape of the PSD, which appears 
not be described by a simple smoothly-varying 
quasi-power-law form. There appears a sharp fea- 
ture at frequency v ~ 7 x 10 _4 Hz which is too 
ill-defined to be claimed as a 'quasi-periodic os- 
cillation', but at least does indicate the presence 
of a sharp transition in the PSD. The significance 
of the potential feature depends on the model fit 
to the PSD. A smoothly-varying quasi-power-law 
would imply a significant excess of power, whereas 
fitting a broken power-law would result in no ap- 
parent excess in power. We do not explicitly con- 
sider the origin of this feature in this paper. As 
previously found, and as commonly observed in 
other AGN, the PSD is higher for the time se- 
ries of higher energy bands at high frequencies, 
although they are more comparable in amplitude 
at low frequencies. 

Lag spectra between various bands have been 
shown for 1H 0707-495 and NGC 4051 bv lFabian etal 
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(|2009h : IZoghbi et alJ (l2010h : lMiller et all (|2010allbh 

and for prev i ous observations of Ark 5 6 4 by 

Arevalo et all (|2006l ): lMcHardv et all (|2007l ): lDe Marco 
(20lJ). The lag spectrum between the 0.4-1 keV 
and 4-7.5 keV bands from the new observation of 
Ark 564 is shown in Fig. [3] with frequency reso- 
lution Alog 10 i/ = 0.15. There are strong pos- 
itive lags up to r ~ 200 s at low freq u encies , 
as previously reported by Arevalo et al.1 (|2006h . 
and negative lags at higher frequencies. We note 
the sharp negative lag of r = —222 ± 78 s at 
v ~ 4 x 10~ 4 Hz accompanied by oscillatory and 
negative lags, r ~ —20 s, at higher frequencies, 
similar to the simple top -hat reverberation t rans- 
fer functions discussed by lMiller et aL ( 2010a ) and 
Miller fc Tureen (|201ll ). Overall, the lag spectrum 
has a closely similar behavior to that previously 
seen in well-studied AGN such as 1H 0707-5495 
( Miller et al.ll2010a ). and that similarity suggests 
that there is a single physical mechanism common 
to all these AGN, a point we return to in section[B] 

The frequency of the sharp negative lag and the 
previously cited mass of Ark 564 are in line with 



et al 



Non-stationarity of the PSD and lag 
spectra 



the fr equency-mass relation of iDe Marco et al 



(|2012j ). The magnitude o f the l ag is much greater 
than the IDe Marco et al.1 (|2012l ) lag- mass relation 
would imply. Howeve r, the hard band used by 



De Marco et al.1 ((2012) is 1-5 keV, very different 
from the 4-7.5 keV hard band used in this analysis. 



2500 



2000 



-500 
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Fig. 3. — Lag spectrum of Ark 564 between the 
0.4-1 keV and 4-7.5 keV bands at frequency reso- 
lution A log 10 v = 0.15. Vertical bars indicate 68% 
confidence region uncertainties, horizontal bars in- 
dicate the range of frequencies incorporated in 
each point. 



In the course of the analysis of the light curve 
peaks described in the following sections, some 
evidence was noticed for non-stationarity of the 
time series (Section I5.6|) . In this section we test 
for non-stationarity in the powerspectrum analysis 
by dividing the set of observations into two halves 
and recalculating the PSD and lag spectra. Fig.@] 
shows the PSD for the soft band 0.4 — 1 keV only, 
and the lag spectra between the soft and hard, 
4 — 7.5 keV, bands, for the first four and last four 
OBSIDs analyzed independently. There is clear 
evidence for the PSD having changed in ampli- 
tude on a timescale of a few weeks between these 
two sets of observations (see Tabled}, and some 
evidence for the shape having changed. 

The lag spectra also appear to have varied: the 
sharp negative feature appears more pronounced 
with a lag r = —306 ± 122 s in the first half but 
less pronounced in the later half, r = — 127±94s: 
this difference on its own is not very significant 
but is worth noting, given the other apparent de- 
partures from stationarity that are seen. The neg- 
ative lags at high frequency in the overall lag spec- 
trum (Fig. [3]) appear not to be present in the ear- 
lier half and to be more oscillatory in the later half. 
We also analyzed each of the eight OBSIDs inde- 
pendently, but given the larger uncertainties did 
not find significant evidence for non-stationarity 
on these shorter timescales. 

Apparent variations in estimated PSDs may in 
principle arise from red noise leak of power below 
the lowest sampled frequency, causing what may 
be t ermed an apparent 'weak non-stationarity' 
(e.g. IVaughan et al.ll2003l ). For the observations 
presented here, the lowest frequencies measured 
are more than two decades lower than the PSD 
break frequency, and the time series variance only 
weakly diverges at the low frequency limit. The 
PSD variations that are observed occur through- 
out the full frequency range, and do not prefer- 
entially occur at low frequencies. It is likely that 
they reflect genuine non-stationarity, rather than 
arising from the presence of inadequately-sampled 
long-timescale variations. This conclusion is rein- 
forced by the changes observed in the time-domain 
transfer function (section f5.6[) . 
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Fig. 4. — (left) PSD of Ark 564 in the 0.4-1 keV band, sampled at frequency resolution A \ogv — 0.05, 
analyzed for the first four OBSIDs (filled circles) and the last four OBSIDs (crosses), (right) Lag spectra 
between the 0.4-1 keV and 4-7.5 keV bands, sampled at frequency resolution Alogf = 0.15, for the same 
data subsets. Vertical bars indicate 68% confidence region uncertainties, horizontal bars indicate the range 
of frequencies incorporated in each point. 



5. Analysis of light curve peaks 
5.1. The goal of the peaks analysis 

While a number of claims have been made for 
observation of X-ray reverberation in AGN, ei 



ther with a measured tim e delay (e.g. iFabian et al 



2009; 
2006; 



Miller et al. l2010bl) or without (jMiller et al. 
Ponti et al.l2012l) . no studies have extracted 



the associated transfer function in an individual 
energy band. For those claims with a measured 
time delay, the lag spectrum has been derived from 
the phases of the cross-spectrum, which itself is 
the Fourier transform of the cross-correlation func- 
tion. Previous studies have measured only the rel- 
ative phases of modes between two bands, which 
does not allow a unique interpretation of the mea- 
sured relative time lags. The three mechanisms for 
time delays proposed to date (Comptonization de- 
lays, reverberation, and fluctuations propagating 
radially across an accretion disk) are all expected 
to be characterized by an energy-dependent trans- 
fer function. To measure the transfer function in 
an individual time series requires a blind deconvo- 
lution algorithm. Extensive work has been done 
on the problem of blind deconvolution, mostly for 
image analysis (see Kundur fc Hatzinakosll996l for 
an overview). However, the presence of significant 
shot noise and gaps in the time series rules out the 
use of many of the simpler techniques. 

Instead of performing a full deconvolution, we 



look at the behavior of the time series in the re- 
gion of flares in activity. The presence of a flare 
brings the emission significantly above the level of 
the complicated 'long-timescale' temporal activ- 
ity. By examining the temporal response to such 
a flare in each energy band, we can obtain an es- 
timate of the shape of the transfer function. In- 
dividual flares suffer too much noise, arising from 
the source variations, to provide useful informa- 
tion on the transfer function, but by averaging 
many flares we can obtain an estimate, as shown 
below. This is similar to the technique of 'flare 
stacking' briefly described, but not significantly 
used, bv lZoghbi et al. I (l2010h . 



5.2. The selection and analysis of light- 
curve flares 

Flares in the time series were selected for the 
analysis and measured by a multi-stage process. 
We considered a model in which the light curve 
comprises a long-timescale 'baseline' variation, 
which we wish to filter out, on which are su- 
perimposed flares, with some unknown temporal 
structure, which we wish to measure. In practice, 
the long-timescale variation itself is likely to be 
composed of many overlapping flares, but we are 
unable to separate those into individual events. 
The more extreme flares are brought above the 
baseline level, and we hope to use this to examine 
their shape. Peaks in the light curve were first 
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selected by median filtering the time series, as 
follows. 

1. The raw soft-band time series was filtered with 
a median filter of length 2R+ 1 time bins, giving 
a first-level baseline and residual time series. 

2. The first- level residual time series was smoothed 
with a smoothing function derived from the soft- 
band PSD. 

3. Peaks of greater than 2a significance were iden- 
tified in the smoothed time series. 

A median filter is better than a linear filter at 
preserving temporal structure in the residual time 
series on timescales less than half the length of the 
filter, but to accurately recover that structure, it 
is necessary that the length of the filter window be 
greater than the size of any features under investi- 
gation. The lengths of the filter for this detection 
stage are given in Table[5J and are significantly 
longer than the FWHM of the detection smooth- 
ing function, as required. 

The choice of smoothing function affects which 
peaks are selected, with a tendency to select peaks 
whose shape and duration are similar to those of 
the smoothing function. To optimize the choice of 
smoothing function, we proceeded as follows. In 
the 'moving average' model of the time series, if 
the underlying time series has a white noise pow- 
erspectrum, the observed PSD is the modulus- 
squared of the Fourier transform of the transfer 
function. Assuming that the transfer function is 
symmetric in the 0.4-1 keV band, we obtained an 
estimate of the shape and length of the transfer 
function by Fourier-transforming the positive, real 
square-root of the PSD. The assumption of sym- 
metry is important as we should avoid introducing 
spurious time lags into the analysis through our 
choice of smoothing function. Because the peaks 
were detected in the residual time series after the 
data had been median filtered, we also applied the 
same procedure to the PSD-derived function to 
obtain the final smoothing function. The FWHM 
of the resulting functions were 500,530,370s for 
each of NGC 4051, 1H 0707-495 and Ark 564, re- 
spectively. While the derived FWHM are rather 
similar between these AGN, this is a reflection of 
the similarity of the shapes of their PSDs, and in 
particular of their break frequencies. 

To test the sensitivity of the results to the 
choice of smoothing function, we also experi- 



mented with a Gaussian smoothing function with 
the same FWHM. We obtained similar results, 
although the detailed shapes of the primary flare 
peaks differed slightly in the two cases, and the 
powerspectra of the selected peaks fsection [5T4^) 
exhibited small differences. 

Having selected a set of peaks in the soft band, 
we then attempted to characterize their tempo- 
ral shapes in the soft and hard bands. To sample 
the long, ks, timescales present in the lag spec- 
tra, a long filter would be required. However, di- 
rectly applying a long median filter to the raw data 
produced unsatisfactory baseline curves, which did 
not follow the light curve on long timescales ow- 
ing to the substantial structure on intermediate 
timescales. The solution adopted was nested fil- 
tering, in effect performing a form of wavelet me- 
dian filtering. The measurement steps continued 
from the peak-detection stage, as follows. 

4. The first-level baseline from the peak-detection 
stage was median filtered with length 4i? + 1 
time bins. 

5. The result was further median filtered with 
length 8-R+l time bins, giving a third-level base- 
line. 

6. This third- level baseline was subtracted from 
the original time series to give a third-level 
residual time series. The three-stage filtering 
process produced baselines that better fit the 
behavior of the light curve than a single stage 
filter of length 8R + 1 time bins. 

7. For each of the peaks detected in step 3, a region 
8R+ 1 time bins in width, centred on each peak, 
was selected from the third-stage residual, to 
measure its shape. 

8. A weighted average of all peak shapes was per- 
formed to give an estimate of the average peak 
shape. 

9. Steps 1, 4-8 were performed on the hard-band 
time series to give an estimate of the average 
hard-band peak shape, using, in step 7, the same 
peak times that were found for the soft band. 

Because we are interested in the mean shape of 
peaks of differing amplitude, each individual peak 
and its neighboring filtered values needed to be 
scaled. This was achieved by dividing each set of 
values, in the soft and hard bands, by the char- 
acteristic height of the soft-band peak. The char- 
acteristic height was found by smoothing the fi- 



8 




2 I ' 1 ' ' ■— I ' " I ' 1 ' I ' Li 

10000 20000 30000 40000 50000 60000 

t/s 



Fig. 5. — A section from the Ark 564 time series. The upper panel shows the soft (0.4-1 keV) band, while 
the lower panel shows the hard (4-7.5 keV) band. Values have been scaled by the overall mean flux in each 
band, and the 68% confidence region due to shot noise is indicated by the grey shading. The solid black 
curve denotes the third-stage median filter baseline. The black dotted line in the upper plot shows the first 
stage median filter baseline used for peak selection. Vertical lines show the times of peaks selected in the 
soft band. 
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nal residual time series with the same smoothing 
function used to select the peaks, to reduce the im- 
pact of shot noise on the central peak value, and 
measuring the maximum value. The final average 
shape in each band was then calculated using a 
weighted average of the individual, scaled but un- 
smoothed, peak shapes. The weights used were 
the total shot-noise variance in the soft band. 

The result is that the hard band peak shape 
is generated using the peak times, scalings and 
weights of the soft band peaks, and is completely 
blind to the hard band light curve. 

5.3. Peak Shapes 

The resulting peak shapes are shown in Fig. [6] 
The peak shapes in each band have been scaled 
such that the height at t = is unity. This al- 
lows for easier visual comparison of the shapes. In 
Ark 564, an approximate comparison of the hard 
band primary peak with the excess gives a unit- 
less area for the primary of 4.6 ± 0.4 over a range 
in time of ±384 s, while the excess has an area 
of 7.6 ± 0.7 over a range in delay time of 640- 
2700 s (the appropriate time range was obtained 
from visual inspection of Fig.[5]and is rather uncer- 
tain). Thus the hard band excess contains about 
1.65 ±0.2 times as many photons as the hard band 
primary peak. In NGC 4051 and 1H 0707-495 
the delayed excess and the primary peak cannot 
be easily separated, making it difficult to evaluate 
the fraction of delayed photons. Visual inspection 
indicates that those fractions are smaller in these 
two AGN compared with Ark 564. These values 
and the shapes obtained are discussed in section[6] 

The parameters used to obtain the peak shapes 
are shown in Table[2] The shapes of the peaks 
were relatively insensitive to changes in the param- 
eters. The time series for AGN 1H 0707-495 and 
NGC 4051 were binned in intervals of 256 s owing 
to the noisiness of the data. Ark 564 has better 
quality data, and was binned at 128 s. Reducing 
the size of the time bins results in noisier peak 
shapes, and difficulties in obtaining clear PSD and 
lag spectra fsection [5T4|) . Increasing the bin size re- 
duces the resolution available on the peak shape. 
The bin sizes were chosen to provide a good com- 
promise between these two factors. The values 
chosen oversample the peak by a factor 2-3, while 
the delayed excess is on a much longer timescale, 
and thus well sampled. 



Since we are only sensitive to timescales shorter 
than the filter window length, the first-stage fil- 
ter length was chosen to be sufficiently long that 
modes at frequencies below the break in the PSD 
were included. The third-stage filter length was 
a factor four larger, such that the low frequency 
modes that show positive lags in the lag spectrum 
were included. In the case of NGC 4051 it was 
found that the filter lengths chosen resulted in sig- 
nificant negative residuals appearing in the final 
mean peak shape, which disappeared when the fil- 
ter lengths were increased by a factor 2, so in this 
case that longer filter length was adopted. 

The main other effect of increasing the filter 
length is to increase the significance of the hard- 
band excess, and to increase the number of peaks 
detected. The latter effect is because as the fil- 
ter length increases, the variations about the base- 
line level increase and more peaks meet the signif- 
icance criteria. The former effect likely arises be- 
cause more low-frequency, positive-lag modes are 
included as the filter length is increased. 

As we were concerned about the possible im- 
pact of overlapping neighboring peaks, we per- 
formed a visual inspection of the peaks located. 
While some peaks do overlap, the overall impact 
is symmetric - for every peak with another closely 
following it, there is a peak with one closely pre- 
ceding. The final impact of the few overlapping 
peaks is minimal. 

5.4. Peak PSD analysis 

We would like to verify that the peak shapes we 
extract are a good estimate of the transfer func- 
tion. To do this, we Fourier transform the peak 
shapes shown in Fig.|6l deduce their power and lag 
spectra and compare those with those from the full 
time series. If the time series is a simple, station- 
ary convolution of an underlying process with a 
transfer function, then we expect the spectra to 
have similar features. 

To find power and lag spectra, the peak shapes 
were padded to a length of 1001 time bins (i.e. 
much longer than the timescale of interest to avoid 
aliasing problems), and Fourier transformed. The 
Fourier modes were then binned to equal sized 
bins in log-frequency. The bin widths were suf- 
ficiently wide to eliminate correlation between ad- 
jacent modes. The low frequency information cre- 
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Fig. 6.— The extracted peak shapes for the three AGN: (top) Ark 564 sampled at 128 s; (centre) 1H 0707-495 
sampled at 256 s; (bottom) NGC 4051 sampled at 256 s. Shaded regions indicate the 68 percent confidence 
region for each plotted point. Shapes are shown for the 0.4-1 keV (dark shading) and 4-7.5 keV (light 
shading) bands. The peak shapes have been normalized so that the amplitude at t = is 1.0. 



Table 2: AGN peak detection and measurement parameters 



AGN 


Ark 564 


1H 0707-495 


NGC 4051 


Sampling time /s 


128 


256 


256 


Peak detection filter length / s 


2688 


2816 


5376 


Peak measurement filter length / s 


10368 


10496 


20736 


Smoothing function FWHM /s 


370 


530 


500 


Peak Significance 


2a 


2a 


2a 


Number of peaks found 


117 


65 


79 
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ated in the padded series below the lowest mea- 
sured values was discarded. The padding was im- 
portant for eliminating aliasing caused by the peri- 
odic boundary condition of the Fourier transform, 
and prevents spurious modes being created within 
the frequency range of interest. 

Bootstrap estimates were made of the statisti- 
cal uncertainties on the Fourier modes. For an av- 
eraged peak shape formed from a weighted average 
of set of n individual peaks, the bootstrap method 
involved generating 1000 realizations, each con- 
sisting of n peaks drawn at random, but with re- 
placement, from the observed set. In a given re- 
alization, some peaks could be present multiple 
times, others not at all. For each realization, peak 
shapes and Fourier transforms were found and the 
variance of the realizations was taken to be the 
variance on the Fourier transform of the actual 
peak shape. The resulting PSDs are shown in 
Fig-El compared with the PSDs of the raw time se- 
ries for the same observations. Because the peaks 
transfer function has been been renormalized, and 
in any case by selecting the highest peaks the am- 
plitude cannot be considered typical of the entire 
time series, the transfer function PSDs have been 
normalized so that the total variance matches that 
of the full time series variance. It can be seen that 
the PSD of the transfer function and of the full 
time series match well for NGC 4051, and also 
match well at low frequencies in 1H 0707-495 and 
Ark 564, but that the transfer function has less 
power at high frequencies in the latter AGN, par- 
ticularly in the soft band. This is discussed further 
in section lGTTl 

Lag spectra between the hard (4 — 7.5 keV) 
and soft (0.4 — 1 keV) bands were derived from 
the transfer functions, and are compared with lag 
spectra from the full time series in Fig.[H Uncer- 
tainties were estimated using the bootstrap pro- 
cess described above. The transfer function lag 
spectra are more noisy, as expected given that only 
a subset of the data has been used, but reproduce 
well the principal features at low frequencies, in- 
cluding the transition to lags around zero at high 
frequency. With the exception of the sharp nega- 
tive lag feature in Ark 564, the transfer functions 
are too noisy to be able to reproduce the negative 
lags at high frequencies. 

Overall, we conclude from the similarity of both 
the PSDs and the lag spectra of the extracted 



peaks and the full time series demonstrate that the 
peaks analysis method does capture both the pow- 
erspectrum and lag information that is obtained 
from a standard Fourier analysis, but now with the 
benefit of being represented as the time-domain 
transfer function in an individual energy band. 

5.5. The dependence on energy band of 
peak selection 

Considering the extreme appearance of the 
hard excess in Ark 564, we undertook a further 
analysis to examine its nature. To do this we also 
investigated time series in a medium energy band 
(1-2 keV). When we perform the peak selection 
in this medium band, the resulting set of peaks is 
different from those selected in the soft band. We 
end up with three sets of peaks - 52 detected only 
in the soft band, 29 detected only in the medium 
band, and 65 detected in both. These peak shapes 
are shown in Fig. [3] The hard-band peak shape 
for this latter category shows a hard excess sim- 
ilar to that found previously. The set of peaks 
found only in the soft band produces an extreme 
hard-excess, with a small primary peak. The set 
of peaks found only in the medium-band shows 
no sign of a hard excess. A close examination of 
the peak positions shows that while no peak was 
detected in the soft-band, the positions are those 
of maxima in the soft-band, just not significant 
enough to make the cut for selection. Whether 
this is due to noise or a systematic weakening of a 
subset of peaks in the soft-band is unknown. This 
analysis is discussed in section l6.2.4l 

5.6. The effects of non-stationarity on the 
peak shape 

In section|4]we have shown that the underlying 
process in Ark 564 is non-stationary. In light of 
this we examine peak shapes obtained from each 
of the two halves of the dataset. The resulting 
peak shapes, shown in Fig.[10l exhibit clear differ- 
ences. Both peak shapes exhibit an excess in the 
hard band at later times than the primary peak. 
The shape from the first half of the data set has 
a larger excess, peaking with an intensity compa- 
rable with that of the primary peak and extend- 
ing over ~ 2ks. A crude comparison of the hard 
band primary peak with the excess gives a unit- 
less area for the primary of 3.9 ± 0.7 over a range 
in time of ±384 s, while the excess has an area 
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Fig. 7. — Power spectra from the extracted peak shapes of the three AGN: (left) Ark 564; (center) 1H 0707- 
495; (right) NGC 4051. The top three panels display the soft (0.4-1 keV) band power spectra, while the 
bottom three panels show the hard (4-7.5 keV) band power spectra. Each panel compares the power spectrum 
of the peak shape we extract (black circles) with the PSD of the full time series (gray squares). Vertical bars 
indicate the 68% confidence region on each point. 
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Fig. 8. — Lag spectra between 0.4 — IkeV and 4 — 7.5 keV bands from the extracted peak shapes (black 
points) compared with that derived from the full time series (gray points) for the three AGN: (left) Ark 564; 
(center) 1H 0707-495; (right) NGC 4051. Vertical bars indicate the 68% confidence region on each point. 
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Fig. 9. — Peak shapes in soft and hard bands for three different selections of peaks in Ark 564: (top) peaks 
detected in both soft (0.4 1 keV) and medium (1 2 keV) bands; (center) peaks detected in only the soft band: 
(lower) peaks detected in only the medium band. Shaded regions indicate the 68 percent confidence region 
for each plotted point. Shapes are shown for the 0.4-1 keV (dark shading) and 4-7.5 keV (light shading) 
bands. 
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of 10.9 ± 1.1 over a range in delay time of 900- 
3000 s. This corresponds to the hard band excess 
containing about 2.8 ±0.6 times as many photons 
as the hard band primary peak. The shape from 
the second half of the data set has a far smaller 
hard excess extending over just ~ Iks. The ap- 
proximate comparison of the hard band primary 
peak with the excess gives a unitless area for the 
primary of 5.5 ± 0.6 while the excess has an area 
of 3.5 ± 0.7 over a range in time of 1500-2700 s, 
corresponding to the hard band excess containing 
approximately 0.6 ± 0.15 as many photons as the 
primary peak (the time ranges were chosen follow- 
ing visual evaluation of the period having the most 
significant delayed excess). 

It is unclear how this effect and that shown in 
section 15.51 are linked. Table [3] lists the number of 
peaks detected in each of the observation periods. 
In the first four observations, peaks detected only 
in the soft (0.4— IkeV) band outnumber those de- 
tected in both soft and medium (1 — 2keV) bands. 
In the last four observations the converse is true. 
Peaks detected only in the medium band appear 
to be uniformly spread. Thus the primary cause 
of the apparent variations in the delayed excess 
is not clear. Does the excess depend primarily 
on the peak spectral hardness ratio, with peaks 
of differing hardness ratios appearing in the time 
series preferentially at differing epochs? Or are 
both the variations in the delayed excess and the 
variations in the hardness of the detected peaks 
separate manifestations of the non-stationary vari- 
ations in the source? 

5.7. The nature of the delayed excess 

The hard excess that is shown in Fig. [6] has been 
obtained by averaging many peaks together. The 
resulting mean shapes for the hard excess could be 
caused by all peaks having a temporally broad, de- 
layed signal with a hard spectrum. Alternatively, 
it might be that an individual primary flare tends 
to be followed by a narrow, secondary peak at a 
random time offset t, where r follows a probability 
distribution such that when averaged, we obtain 
the broad shapes seen. If the latter explanation 
were true, then we would expect that the excess 
in any given time bin would be dominated by a few 
individual peaks. To test for this, we calculated a 
weighted median of the individual peaks for each 
time bin, rather than a weighted mean. If the dis- 



tribution of offset secondary peaks were the expla- 
nation, then we would expect the weighted median 
to fall below the weighted mean. In fact we find 
that the weighted median agrees with the mean 
estimate within the uncertainties, indicating that 
the delayed excess is indeed a smooth continuous 
function that is consistent between all flares. 

6. Discussion 

6.1. Peak shapes and power spectra 

The peak shapes shown in Fig.[B] have several 
common features. The soft-band peak is broadly 
symmetric. The reference time t = is defined by 
this soft-band peak center. The hard-band peaks 
within 256 s of the soft band, and is highly asym- 
metric, exhibiting an excess at later times. The 
nature of this excess varies across AGN, being 
most prominent and extensive in Ark 564. 

The power spectra from the peak shapes 
(Fig. [7]) are in good agreement with those from the 
full time series in NGC 4051 and at low frequen- 
cies in the other AGN, but lack high-frequency 
power in Ark 564 and 1H 0707-495. The break 
frequency in the soft band appears to be the same 
for the peak shape as for the full series, although 
the power at frequencies above the break appears 
to be suppressed in the peak shapes. Three possi- 
ble explanations of this discrepancy are: 

1. There is an uncertainty in the position of the 
peaks detected in the smoothed time series. 
This uncertainty broadens the peak shapes, sup- 
pressing high frequency power. 

2. There might be additional high-frequency vari- 
ations in the full time series, not captured by 
the peak selection. 

3. The existence of the rms-flux relation (e.g. 
McHardvl l2010h may cause aliasing of modes 



such that high-frequency power is boosted with 
respect to that expected in the simple 'moving 
average' model d. 
However, the good agreement between the PSDs 



2 The rms-flux relation may be reproduced mathematically 
by taking a Gaussian time series and multiplying by a si- 
nusoidal function in the time domain. The Fourier trans- 
form of the resulting time series is then expected to be the 
convolution of the Fourier transform of the Gaussian time 
series with the Fourier transform of the multiplying func- 
tion. That convolution aliases power to high frequencies in 
the powerspectrum. 
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Fig. 10. — Peak shapes in soft and hard bands for the first four observations (top) and the last four 
observations (bottom) of Ark 564. Shaded regions indicate the 68 percent confidence region for each plotted 
point. Shapes are shown for the 0.4-1 keV (dark shading) and 4-7.5 keV (light shading) bands. 



Table 3: The number of peaks found in each OBSID, broken down by detection only in the soft (0.4 — 1 keV) 
band, in both soft and medium (1 — 2keV) bands, or only in the medium band. 



OBSID 


0.4 - 1 keV 


both 0.4 - 1 keV 


1 - 2keV 




alone 


& 1 -2keV 


alone 


0670130201 


7 


9 


3 


0670130301 


7 


7 


5 


0670130401 


8 


4 


5 


0670130501 


13 


8 


4 


0670130601 


5 


9 


3 


0670130701 





11 


4 


0670130801 


6 


8 


2 


0670130901 


6 


9 


3 



17 



at the frequencies where the positive lags are seen 
in all three AGN indicates that the peaks analysis 
procedure fully captures the timing information at 
those frequencies. 

When we examine the lag spectra of the trans- 
fer functions we find that the soft band lag is 
consistent with zero (although we note that, had 
we chosen to smooth with an asymmetric smooth- 
ing function, the resulting peaks would have had 
a non-zero lag at some frequencies). Thus, the 
lag on the cross spectrum is primarily due to 
the shape of the hard band transfer function. 
We particularly note the sharp negative lag fea- 
ture at 4 x 10~ 4 Hz in Ark 564 which appears in 
both the Fourier lag spectrum and the PSD of 
the hard-ba n d peak , de monstrating the point o f 
Miller et al.l (|2010al ) and iMiller fe Turner! (|201lh 
that such negative features arise naturally from 
the effects of Fourier analysis of reverberation-like 
transfer functions fsection r6.2.5p . A straightfor- 
ward way to obtain such a negative lag is to have 
a gap in the transfer function between the 'di- 
rect' (zero-lag) and 'rev erberation' (delayed) sig- 
nals ([Miller et al.ll2010al ) , and this seems to be the 
case in Ark 564. However, in other AGN, such a 
gap could be masked by convolution with a broad 
flare temporal response, as appears to be the case 
in 1H0707-495. If the flare temporal response is 
symmetric, the convolution effect does not change 
the cross-spectrum phases, and hence the time lags 
in the lag spectrum, but can reduce the appear- 
ance of a gap in the time domain. 



6.2. Possible explanations of the observed 
temporal behavior 

6.2.1. Intrinsic flare temporal structure 

We now consider possible explanations for the 
lag spectra that have been discussed in the lit- 
erature. At present, we have little idea of the 
mechanism that creates the flares and the gen- 
eral time series variations. It is possible there- 
fore that any temporal structure may be intrin- 
sic to some complex high-energy process and not 
a consequence of propagation of either photons 
or sound through an extended distribution. For 
example, solar flares are generated by a rather 
poorly-understood mechanism, and also show fea- 
tures such as sy stematic spectral va riations in 
hard X-rays (e.g. iKrucker fc Linl 120021) . although 



we note that the solar X-ray spectra and spec- 
tral variations are complex and differ substantially 
from those observed in AGN. Until the primary 
flare mechanism has been definitively established 
in AGN, we would never be able to eliminate this 
possibility. However, we note that the timescale of 
the lags in Ark 564 are substantially longer than 
the duration of the primary flare, and it does seem 
likely that the delayed, hard-spectrum excess seen 
in Ark 564 arises through secondary reprocessing 
rather than through the primary flare process. 

The durations of the primary flares measured 
in section l5~2l of 370 — 500 s allow upper limits to 
be deduced for the sizes of the regions produc- 
ing those flares. The black hole masses are uncer- 
tain (section[T]) , but the flare durations imply light 
crossing-times d < 50 r g . 

6.2.2. Inverse Compton upscattering time delays 

Our analysis disfavors models that explain the 
lag as being due to inverse Compton upscattering 
of photons, in which the hard band variations are 
delayed by the time taken to upscatter photons 
from the soft band. We find no sign of such a de- 
lay of the hard-band peak with respect to the soft- 
band peak in the peak shapes we extract: the posi- 
tive lags appear to be entirely due to the additional 
hard-spectrum delayed component alone. The pri- 
mary flares may be produced by inverse Comp- 
ton upscattering, if the upscattering timescale is 
shorter than the temporal resolution of our analy- 
sis. The hard-band excess has timescales a factor 
> 10 longer, not consistent with arising as inverse 
Compton upscattering delay in the flare genera- 
tion region. 



6.2.3. Positive 
Hons 



lags from propagating fluctua- 



Arevalo fc Uttlevl (|2006l ) have proposed that 



the low-frequency, positive lags in AGN arise from 
fluctuations in the accretion disk that propagate 
radially inwards, with an assumption that a fluctu- 
ation at large radius results in a soft X-ray emis- 
sion spectrum, whereas a fluctuation at smaller 
radius results in a harder X-ray spectrum. This 
assumption is made because the accretion disk at 
small radii is hotter th an at larger radii in the stan - 
dard thin-disk model ( Shakura fc Svunvaevlll973 ). 
The viability of such a model in AGN is open to 
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question. It requires inward-propagating modes 
to dominate the variability of X-rays. Also an 
AGN accretion disk itself is, unlike that of a galac- 
tic black hole candidate, insufficiently hot to emit 
significant hard X-rays. Instead the disk is re- 
quired to act as a source of seed UV emission to 
an inverse Compton upscattering process, and the 
process is thought to operate in the regime where 
the spectral index of the upscattered radiation is 
large ly independent of the s eed photon distribu- 
tion (jTitarchuk fc Huall 19971) . This potential issue 
can be gotten around by allowing radiation from 
different parts of the disk to illuminate different 
Comptonizing regions with different properties. 

This model is also disfavored by the measure- 
ments of Ark 564 presented here: the flaring be- 
havior shows sharp peaks that are temporally co- 
incident at all photon energies, accompanied by a 
delayed, hard-spectrum signal, whereas in the sim- 
plest 'propagating fluctuations' model we would 
expect to see a broad temporal profile in which 
the spectrum hardened with time. Such a model 
would be consistent with the other two AGN, 
whose delayed signals appear as asymmetries in 
the peak shape rather than in the form of a dis- 
tinct component, but if we wish to find a single 
explanation for the lags in all AGN, the observa- 
tion of time-coincident peaks accompanied by a 
distinct, hard, secondary maximum in Ark 564 re- 
mains challenging. 

To explain this, a more sophisticated model of 
'propagating fluctuations' might suppose that a 
flare may be produced at some location, either 
on the accretion disk or in the corona, which is 
directly observed at all energies, but that a dis- 
turbance then propagates into a hard-spectrum 
region such that a delayed, hard signal is seen. 
In Ark 564 there is a distinct gap in the transfer 
function after the primary peak which does not 
have an obvious explanation in this scenario. It 
is worth noting that the gap is responsible for the 
sharp negative lag in the lag spectrum of Ark 564, 
and is a significant feature to be explained (see 
section l6".2.5j) . 

Such a model might, however, be able to ex- 
plain the observation of section l5~5l that the de- 
layed signal is stronger for flares that have soft 
spectra, and absent for those flares that have hard 
spectra. A possible explanation is that the soft 
flares create propagating fluctuations that travel 



into the hard-emitting region, whereas hard flares 
might be created directly in that region. 

6.2.4-. X-ray reverberation 

The final model co nsidered here is the rever- 
beration explana t ion of Miller et al. (2010a. b) and 
Miller fc Turner! (|201lf) . In those papers, it was 
shown how a long-timescale delay, such as might 
be produced by reverberation from gas a few tens 
to thousands of gravitational radii from the X- 
ray source, could naturally produce low-frequency 
positive lags, and also potentially generate high- 
frequency negative lags as a result of the shape of 
the transfer functions in the two b ands being cross- 



correlated ([Miller et all l2010at iMiller fc Turner 



20111 ). Such models are consistent with the anal- 
ysis presented here, and in particular with the 
delayed, hard excesses that are seen in all three 
AGN, but particularly prominently in Ark 564. 
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Fig. 11. — Example model spectrum from the 
uniform spherical scattering model described in 
the text. The upper black curves shows the total 
spectrum received by the observer. The directly 
received component (red curve) and delayed, scat- 
tered component (blue curve) are also shown. 

The fraction of hard-band photons in the de- 
layed excess to photons in the primary peak is 
very high, with a mean value of 1.65 ± 0.2, but 
apparently being as large as 2.8 ± 0.6 in the first 
half of the observations f section l5~6"| . Within the 
reverberation hypothesis, this could arise if the ob- 
server's sight-line to the primary flare site were 
through material that is moderately optically- 
thick to Compton scattering, so that most photons 
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from the flare are scattered and appear in the de- 
layed excess. The dependence on flare hardness 
could arise naturally in ionized gas as, for a wide 
range of values of the ionization parameter, the 
photoelectric opacity in the 1 — 2 keV band can be 
larger than that in the 0.4—1 keV band. Thus in 
the first half of the observations, we could be look- 
ing through a higher optical depth of scattering, 
absorbing gas than in the later half, leading both 
to softer flares around 0.4 — 2 keV and to a larger 
delayed excess in the 4 — 7.5 keV band. 

As a 'toy model' example of this, we have cal- 
culated the spectrum expected from a uniform 
sphere of gas near the surface of which is embed- 
ded a point source of X-rays with a powerlaw con- 
tinuum. Spectra have been calculated by a 3D 
Monte-Carlo radiative transfer algorithm with full 
treatment of Compton scattering by cold, free elec- 
trons, but assuming that photoelectric opacity is 
provided by gas with an invariant ionization pa- 
rameter j_wdiCTe_o_rja£^esliaye been obtained from 
XSTAR (jKallman et al.ll2004l) . No thermal or ion- 
ization balance calculations have been made in- 
side the simulation and line and radiative recom- 
bination continuum emission has been ignored. 
The effects of thermal and macroscopic Doppler 
shifts have been approximated by smoothing with 
a Gaussian velocity dispersion a — 0.05c, appro- 
priate for material at radii r ~ 100 r g . This simple 
model is not expected to be able to fit data in any 
detail (and we stress that this does not represent a 
spectral fit to any of the sources presented in this 
paper), but it does illustrate the principal features 
of having a high fraction of scattered light with a 
hard spectrum. Fig. 1111 shows an example for ion- 
ization parameter log 10 £ = 2.1, electron column 
density N c — 4 x 10 23 cm -2 and mean absorber 
hydrogen column density TVh — 4 x 10 22 cm -2 in 
front of the X-ray source (in practice achievable 
by having a two-phase medium of hot, free elec- 
trons within within are embedded denser clumps 
of lower ionization gas) . The source has a power- 
law index V — 2.5. The depth of the source within 
the sphere is 1/50 the sphere radius. This con- 
figuration gives a ratio of scattered to direct light 
within the 4 — 7.5 keV band of 1.65 as in Ark 564. 

An alternative explanation for the high ampli- 
tude of the delayed excess might be that the re- 
verberating region is itself sufficiently hot that ad- 
ditional photons are inverse Compton upscattered 



into the hard energy band by the reverberating 
gas. Such models may also naturally produce a 
high fraction of X-rays in the delayed excess, as 
there is a large source of soft photons that may 
be suitably upscattered. We would expect the re- 
sulting spectra to be largely featureless, generally 
consistent with the observations of Ark 564. Mod- 
els that seek to link the transfer function behavior 
and the spectral properties of Ark 564 will be dis- 
cussed by Giustini et al. (in preparation). 

The scale-free example model presented above 
may be scaled in size to produce the observed time 
lags. However, we do not expect that a simple 
uniform sphere will reproduce the detailed shape 
of the transfer function, for which more extensive 
radiative transfer modeling of clumpy, anisotropic 
distributions will be required. 

The differing spectra of flares inferred in section 
15.51 need not imply rapid variations in ionization 
parameter or optical depth. Instead, it could be 
that flares are emitted from multiple regions, each 
with their own value of optical depth. The longer 
timescale changes described in sections l4l fe 15 .61 are 
of order weeks (the time between observations). 

6.2.5. Negative lags 

Negative lags are seen in the lag spectra of 
1H 0707-495 and Ark 564, but not in NGC 4051. 
In Ark 564 there is a sharp negative lag which is di- 
rectly linked to the gap in the transfer function be- 
tween the primary peak and the delay ed signal (see 



the dis cuss ion of transfer function s bv lMiller et al 



2010a! and [Miller fc Turnerll201lh . T his seems to 



confirm one of the basic points of iMiller et al 



(|2010al ). that sharp negative oscillations may be 
caused simply by the shape of the hard band trans- 
fer function, and not by any requirement for a sec- 
ondary scattering process in the soft band. 

However, the peaks analysis presented here 
lacks sufficient signal-to-noise to detect the oscil- 
latory, negative lags found at higher frequencies 
in the lag spectra of these two AGN. Persistent 
negative lags that extend over a wide range of 
frequency are difficult to explain if it is only the 
har d band that has an extended transfer func tion 
(see IMiller et al.ll2010ai IMiller fc Turnerll201ll and 
Zoghbi et al.l 201 lh and it seems that an extended 
transf er function may be needed also in the soft 
band ((Miller et al.ll2010al) . In the case of Ark 564, 
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those high-frequency negative lags are low ampli- 
tude, oscillatory and may also be non-stationary 
(section|4j . If it is indeed required, our analy- 
sis cannot detect such a soft-band transfer func- 
tion, as it is limited by the attainable signal-to- 
noise, but we should nonetheless discuss the re- 
quirements that the possible models would need 
to satisfy in each case. 

The 'propagating fluctuations' model does not 
naturally explain those negati ve lag s in 1H 0707 



495, and iFabian et al.l (|2009h and IZoghbi et al" 



(|201(j l201lh tiave proposed that these arise as 
a separate process, in which a compact X-ray 
source, close to both the black hole and the ac- 
cretion disk, illuminates the disk and produces 
X-ray emission. With a suitably-tuned reflection 
spectrum, enhanced by gravitational bending of 
light near the black hole, it was proposed that 
this may result in a delayed reverberation sig- 
nal in the soft band. Some of t he difficulties of 
this m odel have been discussed bv lMiller fc Turner 
(|2011[ ) and are not discussed further here, apart 
from noting that the model requires (i) a very 
high abundance of iron ( at least a factor 1 1 super- 
solar, |Zog^reEal][20i3), (ii) the X -ray source to 
be located within a few gravitational radii of the 
event ho rizon of a maximally-rotati ng black hole 
(see also IChainakun fc Young! 120121 ) and (hi) the 
AGN to be significantly super- Eddingtorfl Thus 
an inverse Compton upscattering source, required 
to produce the powerlaw hard X-ray spectrum and 



' P' . 
postulated by lArevalo et al. (2006), would need 



to be extremely close to the black hole and to 
produce short-timescale flares with long-timescale 
variable emission that is both delayed and signifi- 
cantly harder than the prompt flare spectrum. 

As Ark 564 also shows similar negative lags 
to 1H 0707-495, the complex spectral model in- 
voked for 1H 0707-495 should also be required 
for Ark 564, yet the 0.4-10 keV spectrum of 



3 The blac k hole mass for 1H 0707-495 of 2 X 10 6 M Q as- 
sumed by [Zoghbi et al. I d2010t) led to those authors asso- 
ciating a 30 s negative lag with reflection from a radius 
2 r . Given the inferred bolometric luminosity l lLeighlvl 
l2004h such a black hole would be about a factor 5 super- 
Eddington. If the black hole mass were higher, the radius of 
the reflecting region would be correspondingly smaller, and 
becomes implausibly close to the event horizon for black 
hole masses > IO^Mq. We note that the black hole mass 
of 1H 0707—495 has not yet been determined by optical 
reverberation mapping, at the time of writing. 



Ark 564 appears to differ substantially from that 
of 1H 0707-495 (see Giustini et al., in prepara- 
tion) . 

iMiller et al.l (|2010aj) proposed that, in the re- 
verberation model, some soft-band X-rays scat- 
tered with short delays in a continuous medium 
could reach the observer, but that long-delay scat- 
tered soft X-rays would be absorbed within the 
scattering medium owing to the longer path-length 
traversed, leading to a soft-band transfer func- 
tion that was extended over short timescales only. 
While such a scenario is possible, photoelectric ab- 
sorption can only remove photons from an X-ray 
energy band, implying that the soft-band transfer 
function should always be lower in amplitude than 
the hard band transfer function, so this mecha- 
nism might require fine-tuning to obtain an ap- 
propriate cross-correlation signal. 

If we follo w Fabian et al. ( 20091 ) and Zoghbi et al 
( 2010l 2011 ) , and suppose that there are two sep- 
arate mechanisms that operate at low and high 
frequencies, an alternative possibility is that the 
short timescale (< 30 s) negative lags are actu- 
ally a feature of the emitting region itself. Lags 
on that timescale are significantly shorter than 
the duration of the flares (~ 400 — 500s), and so 
are more plausibly associated with the primary 
process than are the ks-timescale delayed signals 
seen at low frequencies. If the flare-emitting re- 
gion were slightly larger at soft energies than at 
hard energies, soft-band flares would have PSDs 
that cut off at lower frequency than hard-band 
flares (or, equivalently, slightly broader peaks in 
the peaks analysis), as observed, and would have 
variations delayed by time At ~ Ar/c, where 
Ar is the difference in apparent size of the flare 
region between the soft and hard X-ray bands. 
Such a size differential could arise either from the 
effects of photoelectric opacity or from the energy- 
dependent diffusion of electrons from a flare ac- 
celeration region, given that radiative lifetimes 
due to inverse Compton scattering increase with 
decreasing electron energy. 

7. Conclusions 

This paper has explored the origin of the PSDs 
and time delays in three, highly X-ray variable, 
narrow-line SeyfertI AGN: NGC 4051, 1H 0707- 
495 and Ark 564, using previously-published data 
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for the first two and a new long XMM-Newton ob- 
servation of Ark 564. The lag spectrum of Ark 564 
shows strong positive lags (hard photons lagging 
soft photons) for Fourier modes wit h low frequen- 



cies, consistent with previo us work (jArevalo et al 



20061: iMcHardv et al.ll2007l ). It also shows a sharp 
negative lag at 4 x Hz and oscillatory, slightly 
negative lags at higher frequencies. The PSD of 
Ark 564 has a complex shape that includes non- 
smooth structure, but which cannot be clearly la- 
beled as a quasi-periodic oscillation. We have also 
found evidence for non-stationarity of the PSD 
and lag spectra of Ark 564 on timescales of weeks. 

To better understand the lag spectra of these 
AGN, we have developed a method for estimating 
the transfer function of AGN X-ray time series in 
individual bands of photon energy, which allows us 
to investigate physical models of time lags with- 
out the limitation of using cross-correlation be- 
tween bands. By modeling the AGN variations as 
a 'moving average' time series process, the method 
obtains the mean shape of emission peaks in the 
light curve, after a nested median filtering pro- 
cess to remove long-timescale variations has been 
applied. All three AGN show an excess of de- 
layed emission in hard X-rays, but not soft X- 
rays, on ks timescales following an emission flare. 
The delayed emission detected in the peaks anal- 
ysis is able to reproduce the lags between hard 
and soft photons at low frequencies that are de- 
duced from Fourier analysis of the time series. The 
sharp negative lag feature in Ark 564 is seen to 
arise from the shape of the transfer function in 
the hard band, and not from time delays in soft 
X-rays, as discussed in the context of mod els of 
X-ray reverberation by Miller et al. ( 2010bh and 
Miller fc Turner] (|201lh - 

In Ark 564 the delayed hard excess is partic- 
ularly prominent and appears to be a distinct 
secondary reprocessing of the flare photons, as 
expected from reverberation caused by Compton 
scattering of photons by circumnuclear material. 
A high optical depth to Compton scattering is im- 
plied by the high fraction of photons in the delayed 
excess. The amplitude of the excess appears to 
vary both with time and with the hardness of the 
selected peaks, implying that we are observing a 
complex, rapidly-varying region of this AGN. The 
temporal coincidence in all three AGN of the pri- 
mary peaks of emission between hard and soft pho- 



tons is inconsistent with the time lags being caused 
by inverse Compton upscattering delays, and the 
shape of the Ark 564 excess is difficult to recon- 
cile with a model in which time lags are caused 
by fluctuations propagating across an accretion 
disk. Given the estimated black hole masses in 
these AGN, the transfer functions we estimate are 
broadly consistent with reverberation from mate- 
rial 10s to 100s of gravitational radii from the cen- 
tral source and provide further evidence for the ex- 
istence of a substantial covering fraction of scatter- 
ing material in the central regions of type I AGN. 
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